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We consider the constraints on string networks with junctions in which the strings may aU be 
different, as may be found for example in a network of {p,q) cosmic superstrings. We concentrate 
on three aspects of junction dynamics. First we consider the propagation of small amplitude waves 
across a static three-string junction. Then, generalizing our earlier work, we determine the kinematic 
constraints on two colliding strings with different tensions. As before, the important conclusion 
is that strings do not always reconnect with a third string; they can pass straight through one 
another (or in the case of non-abelian strings become stuck in an X configuration), the constraint 
depending on the angle at which the strings meet, on their relative velocity, and on the ratios of 
the string tensions. For example, if the two colliding strings have equal tensions, then for ultra- 
relativistic initial velocities they pass through one another. However, if their tensions are sufficiently 
different they can reconnect. Finally, we consider the global properties of junctions and strings in 
a network. Assuming that, in a network, the incoming waves at a junction are independently 
randomly distributed, we determine the r.m.s. velocities of strings and calculate the average speed 
at which a junction moves along each of the three strings from which it is formed. Our findings 
suggest that junction dynamics may be such as to preferentially remove the heavy strings from the 
network leaving a network of predominantly light strings. Furthermore the r.m.s. velocity of strings 
in a network with junctions is smaller than l/\/2, the result for conventional Nambu-Goto strings 
without junctions in Minkowski spacetime. 



I. INTRODUCTION 

The evolution of cosmic string networks and their 
observational consequences are attracting a great deal 
of interest, particularly since they may lead to indi- 
rect observational tests of superstring theory. The an- 
nihilation of two branes at the end of brane inflation 
[U, 0, H, 0, [Hi is thought to lead to the formation of cosmic 
superstrings which can be cither fundamental F-strings, 
Dirichlet D-strings, or (p, (/j-strings, bound states of the 
two d, 0, H, 01 ■ Most importantly, the predicted tensions 
of these strings are not only compatible with current ob- 
servational bounds (G/.i<2.3 x 10"'' using the third year 
WMAP data but they also lie in a window which 

may be testable with the future LISA gravitational-wave 
detector [lil, [H, [ll . 

Whilst cosmic superstrings share many of the prop- 
erties of standard grand unified (GUT) cosmic strings 
[3, [3 J they differ in some important respects. First, 
when two F-strings (or two D-strings) intersect, they do 
not necessarily 'intercommute' — or exchange partners 
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— with probability P = 1 [T^, ■ (An exception for lo- 
cal cosmic strings is when they intercommute at speeds 
very close to the speed of light [ll|.) The effect of reduc- 
ing P is to increase the density of strings in the final scal- 
ing solution and hence the gravity wave signal, though 
the exact manner in which this happens is under debate 
[ni,[i3,[il,[ll,[13l- Second, the different kinds of strings 
in a cosmic superstring network can meet at junctions 
[2H . Thus when an F and D-sting intersect they cannot 
intercommute, but rather two junctions are formed and 
the original strings become connected by a third {p, q) 
bound state. Junctions also exist in non-abelian string 
networks, for which the fundamental group Tri{M) of the 
vacuum manifold A4 is non-abelian. 

Even though several authors have addressed both an- 
alytically and numerically the cosmological evolution of 
string networks with junctions [H, HI [H, [H, [H, Hi, 
[27l . [28j , the late-time behavior of the network is still not 
fully understood, particularly when there are different 
types of string. One possible outcome is that the junc- 
tions play a minor role and the network reaches a scaling 
solution similar to that of a GUT string network. In 
that case it is possible to make predictions for cosmic 
superstring networks by extracting the information from 
the GUT string simulations, with suitable allowance for 
P ^ 1 and in principle for different tensions. A sec- 
ond possibility is that the presence of junctions makes 
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the network freeze and eventually dominate the energy 
density of the universe. Similar questions have been 
addressed for networks of domain walls with junctions 
[2§,[3S[3l|,[33,[3l,[3|. 

In a previous letter [s^ we studied the dynamics of 
junctions in a local string network in which the indi- 
vidual strings have no long-range interactions and are 
well described by the Nambu-Goto action. We set up 
the equations of motion for three strings of tensions i^ij 
{j = 1,2,3) meeting at a junction at position X(t), and 
were able to solve for the dynamics of X(i) as well as 
to determine how the junction moved along each of the 
strings. A related approach has also been developed in 
the context of representations of baryons as p ieces of open 
string connected at one common point [36l . [sj ■ Having 
constructed this formalism, we initially presented a sim- 
ple highly symmetric exact oscillating loop solution. As 
in the case of ordinary cosmic strings, the existence of 
exact loop solutions may be important in analyzing the 
likely behavior of loops in general. Loops of strings with 
junctions generally evolve rather differently to standard 
string loops; they do not oscillate periodically and initial 
analysis indicates that the number of cusps on the loops 
may be rather different than for conventional Nambu- 
Goto loops. We intend to return to these points in a 
future publication. 

Our second result concerned the intersection of two 
strings with equal tensions /ii, meeting at an angle a, 
but with equal and opposite velocities v. When these 
strings collide, for them to exchange partners and become 
linked by a third string of tension /is requires a strong 
kincmatical constraint. We also discussed the case of 
non-abelian strings, which may become joined by a new 
string lying along the direction of motion of the two initial 
strings, and we found the kinematic constraint for this 
case too. 

We expect that these kinematical constraints could 
have significant consequences on the evolution of string 
networks with junctions: if the relative velocity of the 
strings is large, no junction will form and the strings will 
pass through each other. This in turn means that fewer 
loops will be formed, and hence that the network would 
radiate less energy. It may thus be important to include 
these constraints in analytic and numerical models of net- 
work evolution with junctions. 

In this paper we extend the work of [s^ . and focus on 
three aspects of junction dynamics which we believe will 
be important to understand the global properties of the 
network. After a review of the equations of motion de- 
scribing junction evolution in section|TTl we first consider 
the propagation of small amplitude waves across a static 
junction formed of three strings with generally different 
tensions fsection llll)) . Such small waves will inevitably be 
present on strings in a network, and will radiate gravita- 
tionally. We determine the fraction of energy transmitted 
and reflected across the junction. Then, in sections IIVI 



andlVl we generalize the kinematic constraints of 35[ to 
the case in which the two colliding strings have different 



tensions fii and ^2 ■ Since the initial velocity of strings is 
crucial to determining the importance of the kinematical 
constraints, in sections IVll and IVlIl we consider the global 
properties of junctions and strings in a network. Assum- 
ing that, in a network, the incoming waves at a junction 
arc randomly distributed, we determine the r.m.s. veloci- 
ties of strings and calculate the average speed at which a 
junction moves along each of the three strings from which 
it is formed. Our findings suggest that junction dynamics 
may be such as to preferentially remove the heavy strings 
from the network. Finally we conclude in section IVIIII 



II. REVIEW OF THE BASICS 

We begin by briefly reviewing the equations of motion 
for strings that form Y junctions obtained in [3^. 

We choose the temporal world-sheet coordinate to be 
the global time, t — t, and use the conformal gauge in 
which the spatial coordinates x((t, t) satisfy 



X • x' = 0, 



x2+x'2 



1, 



(1) 



where x = 9tx and x' = d^^. The action for three strings 
of tensions fij (j = 1, 2, 3) meeting at a junction is 

j 

+ dtf,(t)-[x,(,s,(t),t)-X(t)] (2) 

j 

where X(t) is the position of the vertex, are Lagrange 
multipliers, and the Sj(t) are the values of the spatial 
world-sheet coordinate at the vertex. It is straightfor- 
ward to generalize this action to include other vertices, 
present for example in loop configurations. 

Varying Xj yields the wave equation, with solution 



X, (^T, t) = - [a, {a + t) + b, (a-t)], (3) 



where, in order to satisfy the gauge conditions ([T]), 



a-=b- = L 



(4) 



The Lagrange multipliers impose boundary conditions 
that may be written 



ajisj +t) + hj{sj -t)^ 2X(i), 



(5) 



while varying X gives a relation between the Lagrange 
multipliers which can be reduced to 



^Mj[(l + s,)a^ + (l-s,)b;]=0. 



(6) 



Eliminating the a^- using the time derivative of ^ then 
yields 



(7) 
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where 



where 



11 = III + /i2 + 



(8) 



At the junction, the amphtudes b^- of the incoming 
waves are known from the initial conditions; those of the 
outgoing waves, a' , together with the position of the ver- 
tex, are given by |5]) and ([6]), provided the Sj are known. 
Finally, the latter may be found by solving the constraint 
equations a^-^ = 1. The result depends on the string 
tensions and on the angles between the unit vectors b^-. 
Relative to [s^, we introduce a slight change of notation, 
writing for example ci = b2 • bg (denoted by C23 in [1^). 
It will be useful to introduce the quantities Vj defined, 
for example, by 



t^i = M2 + M3 - Ml- 
Then the value of s^. is given by 

MA/fe(l-Cfc) 



(9) 



(10) 



where the Alj are defined by Afi = along with two 
similar equations, and 



A^ = ^A/,(l-c,). 



(11) 



Notice that since the Sk < 1, it follows that the string 
tensions satisfy triangle inequalities, 



V, > 0. 



(12) 



Finally, it follows from pO| that the Sj satisfy the energy 
conservation condition 



(13) 



III. SMALL-AMPLITUDE WAVES 

A first interesting application of these results is to the 
reflection and transmission of small-amplitude waves at a 
junction. Small amplitude waves are expected to build up 
on strings as a result of self-intersections and interactions 
with other strings, and they will radiate gravitationally. 
It is important to understand how they propagate across 
a junction. 

Consider three straight, static strings, for which 



a'j = h'j = {cos9j,sm6j,0) = e 



J' 



(14) 



with the junction at X = 0. The equilibrium condition 



determines the angles between them through pOl 



for example, we have 

ci = cos(6l3 - 62) = 
sin(03-^2) = 



? 2 

MI - MS - 

2M2M3 



2M2M3 



(15) 



(16) 



Note that as a consequence of p2)) . A is real. 

Now suppose there is a small incoming perturbation 
on the first string, so that b'^ = ei + ^b^, with 



Sh[{s) = efi cos ks, 



(17) 



where e is a small dimcnsionless parameter and fi is a 
unit vector satisfying ei • fi = 0; we set Sh'2 = (Jbg = 0. 

The simplest case is where fi = (0, 0, 1) since then 
all Scj = 0, and therefore 6sj = so that the junction 
remains at the same value of a on all strings. In that case, 
the outgoing waves are also in the z-direction. From ([7]) 
we find 



and hence 



(5X(t) = — — eficosfct 
M 



Sa[{s) = — eficosfcs, 
M 

(5a' (s) = ,5a'(s) = -^efiCosfcs. (19) 
M 

It is then straightforward to find the fraction of energy 
transmitted along the strings of tension ^^2^3 and reflected 
along the string of tension fii. Since 



(18) 



we find 



Ri 



5E, = ^ 
J 4 



T2 



dcr{5a.f + Shf), 



, MiM2 
2 ■ 

M 



1^3 



,MiM3 
M 



(20) 



(21) 



For instance, suppose that Mi ^ M2 ~ M3i so that the 
initial perturbation is along the light string. Then, as 
expected, almost all the energy is essentially reflected off 
the junction. 

The results are more interesting if we take fi in the 
xj/-plane, because then the values of Sj do oscillate and 
the junction no longer stays at a fixed position on each 
string. Let us set 



fj = ( — sin 9j , cos 6'j , 0) , 



(22) 



so that Gj • fj = for each j, and again take 5h'^ of the 
form (fTT]) . with Sh^ = 5h'r^ = 0. Then we have 



Sci = 0, 

SC2 

<5c3 



2miM3 
A 



-e cos kt, 



(23) 



2miM2 



-e cos kt. 



It follows that 



Ssi = — .^'^•^ ^ ecosfct. 



6S2 



A 

A 



e cos kt = —Sss- 



(24) 
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Despite this difference, the expressions for the outgoing- 
wave amphtudes are very similar, except that they are 
no longer all in the same direction. We find 

SXit) = e cos kt i^itl^£^e, ^^fi) (25) 

from which it follows that 

Sa'i{s) = -efi cos ks 

Sa'^{s) ^ ^ei2 cos ks, (26) 
M 

Sa'^is) — — ^efscos/cs. 
M 

The fact that the amplitudes of the outgoing waves are 
the same in magnitude independent of the orientation of 
fi is at first sight remarkable, but actually it follows from 
the fact that we can regard the waves as representing 
massless particles propagating along the strings. The 
amplitudes of the outgoing waves can be derived from 
conservation of energy and momentum. It follows that 
the reflection and transmission coefficients are as given 
in (PT|) , a result that has also been obtained in the more 
general setting of multivolume junctions by [3^. 

Comparing ([M]) with one interesting point 

emerges: in the first case, sl'i is in phase with h'l, whereas 
in the second case it is in anti-phase. This means that 
if the incoming wave is at an intermediate angle, the re- 
flected wave is tilted in the opposite direction. 

IV. COLLISIONS OF STRAIGHT STRINGS 

As summarized in the introduction, in [s^ we dis- 
cussed the collision of two straight strings of equal ten- 
sion, /ii = fi2, and derived kinematical constraints on 
such a process. Here we wish to extend this discussion to 
the case of unequal tension, and hence it will be useful 
to define 

^i+ = Hi+H2, ^1-/12. (27) 

Since > 0, /X3 is bounded by 

M-<M3<M+- (28) 

This extension is non-trivial as the problem now lacks the 
symmetry associated with the equal tension case, and 
as a first step one must determine the orientation and 
velocity of the joining string after the collision. This is 
the aim of the present section. 

Consider two strings of tension /ii and ^2 parallel 
to the xy-plane but at angles ±a to the cc-axis, and 
approaching each other with velocities ztw in the z- 
direction. Before the collision (for < < 0, see Fig. [T]) 
we take 

xi,2(o',0 = (— 7^"^crcosa,=F7~"^crsina,±ut), (29) 




FIG. 1: Two colliding strings, of unequal tension, joined by a 
third string (an x-\mk). 

where 7"^ = Vl — i'^ . Thus 

a']^ 2 = (—7^"'" cos a, ^7""'" sin a, ±w), 

bj 2 = (— 7^"^ cosa, =F7~"^ sina, ^F'^)- (30) 

At the collision, i = 0, we suppose that the strings inter- 
change partners. Now, since 7^ 0, we can no longer 
conclude that the third - joining - string formed after 
the two original strings exchange partners, will be static 
and that it will lie either along the x- or y-axis. It is 
still true that this third string must be parallel to the 
xy-plane and that there are two orientations for it, de- 
pending on which ends of the original strings are joined 
to each other. To be specific, let us assume that the two 
ends in the positive- a; region are joined at a new vertex, 
X, while those in the negative- a; region are joined to the 
other end of the new string which has tension jij,. We 
call this an x-link — see Figure [T] (The corresponding 
y-link is obtained by replacing a — > (7r/2 — a) and will 
be discussed below.) 

The new string of tension /i3 thus lies at an angle 9 
to the X axis, and moves in the z-direction with some 
velocity u; thus 

X3(CT,i) = (7~^o'cos6',7~^(Tsin0,ut), (31) 

where = 1/v 1 — u'^- We now determine and u. It 
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follows from (|3T|) that 

b^ = (7-1 cos 0,7-1 sin (32) 

and thus 

ci = — 7-^7--^ cos(a + 6') — uu, 

C2 = — 7j7^7~^ cos(q; — 6') + uw, (33) 

C3 = 27-^ cos^ a - 1. 

The position of the vertex is of course X(<) = X3(s3(t), t), 

I 



so from ((3T|) . 

X= (S37-' cos 0,537^1 sin 0,7.). (34) 

Thus, from ([7]) and pO]). we obtain 

MX + Ah{l - C3)H = -Mi(l - ci)b'i - il/2(l - C2)b^, 

(35) 

whose three components are 



[7Ws3 + ^/3(l-C3)]7,7^cos6i = [Afi(l-ci) + A/2(l-C2)]7"^cosa, 

[Als3 + M3(l-C3)]7;7'sin0 = [A/i(l-ci)-Af2(l-C2)]7~'sina, (36) 
[M - Ah{l - c:i)]u = [A/i(l-ci)-Af2(l-C2)]7;. 



Dividing the y-component of ([36]) by the x-component, 
and comparing with the ^-component gives 



tan^ _u _ Mi{\ - ci) - A/2(l - ca) 
tana ^ v ^ Mi{\ - ci) + A/2(l - C2) ' 



(37) 



The first equality here allows us to eliminate 0. The 
second, combined with ([55]) . then allows us to obtain an 
equation for w, which simplifies to a quadratic for (see 
Appendix A): 

/i?_(sin^ q;)m^ + [^3(1 - i;^) + ^^(w^ cos^ a - sin^ q;)]u^ 



22 2 
/ill) cos a 



0. 



(38) 



This equation always has one positive root which then 
determines 9 from ([57]) . Notice that when /i_ = as in 
[3^ . then M = and = 0: this is the symmetric case. 
On the other hand, in the limit when the bound in pS]) 
is saturated, i.e., — > /ig, then u — > u and 9 a. Also, 
when ^ 1 for /x_ > 0, then u ~* 1 and again 9 ^ a. 

The intercommutation produces kinks on the original 
strings, moving along them at the speed of light; they are 
in fact at the same positions as in the equal tension case, 
namely at 



^1.2 = (7 1 cos a, ±7 ^sina,iv)t. 



V. KINEMATIC CONSTRAINTS 



(39) 



Having found the orientation and velocity of the con- 
necting string, the physical condition that the junctions 
must move apart imposes the requirement 



S3 > 



(40) 



where S3 is given in pi7|) . This in turn implies kine- 
matic constraints on the values of v and a of the colliding 



strings for which an x-link can be formed. Replacing a 
by 7r/2 — a gives a similar condition for the formation of 
a y-link. 

When fi- — 0, constraint (HDl) reproduces the results 
of [11], namely 



/ /i37 \ 

a < arccos - — ) (a;— link) 
/ ^37 \ 

a > arcsin (j/— link). 



(41) 
(42) 



If these bounds are not satisfied, then abclian strings 
presumably pass through one another. Essentially this 
same result, in the special case of all tensions equal 
/ii = //2 = M3i obtained earlier in the context of 
hadronic strings [39j . and also when /i_ = from a dif- 
ferent starting point in the context of Type-I strings [i^ . 

For non-abclian strings, as discussed in [s^ . there is 
a third possibility namely that a link forms in the direc- 
tion of motion of the two initial strings (a z-link) . Indeed 
if the string fluxes on strings 1 and 2 do not commute, 
then for topological reasons the strings cannot simply 
pass through one another: they may only do so with the 
formation of a third string. Note that in this case we 
do not have strings with three different tensions meeting 
at a vertex: at one end of string 3 it is attached to two 
segments of string 1, and at the other to two segments 
of string 2. Hence the linking string in this case does lie 
along the z-axis, simplifying the analysis. In effect, we 
have two vertices of the type discussed in [3^ . The tri- 
angle inequalities now require that both 123 < 2ni 
and /i3 < 2/^2- Furthermore the kinematic constraint in 
this case is that the total length of string 3 must increase; 
we must add the rates from the two ends. It is not abso- 
lutely necessary that S3 at one particular vertex should 
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FIG. 2: Kinematic constraints for ^_ =0. Allowed regions 
for a;-links are to the left of the full curves; for y-links to the 
right of the dashed curves; and, for z-links in the non-abelian 
case, above the horizontal dotted lines. The values of ^3 are 
1.4 (red), 1.2 (blue), 1.0 (black). Allowed regions are shaded 
for the /i3 = 1.4 case. 



be positive. The required condition is 



2^.1 - fi3V 2fi2 
Solving for v gives the limit 



M3« 



> 0. 



V > 



2(Atl + /X2)M3 



(43) 



(44) 



These different constraints are summarized in Figures 
[5] and [3] showing the {a,v) plane for /i = 3. Fig. [2] shows 
the three constraints for the case /i_ = (equations (|¥T|) . 
(|42|) and p4l) ) and various values of fi^ (distinguished by 
color). For non-abelian strings there are three different 
regimes: i) if /is > v^/ii there is a kincmatically forbid- 
den region for all a; ii) for 2/ii/\/3 > /is > \/2/ii only a 
restricted range of a are forbidden; in) for /ia < 2/ti/\/3 
the whole {a,v) plane is included in at least one of the 
allowed regions. Allowed regions for the formation of an 
x-link are to the left of the solid lines, those for a j/-link 
to the right of the dashed lines, while for a z-link in the 
non-abelian case they are above the horizontal dotted 
lines. 

Fig. [3] shows the allowed regions for a number of cases 
with > 0. (The sign of does not affect the limits.) 
The allowed x-link regions are to the left of the curves, 
and those for a z-link above the horizontal lines. (The y- 
link regions can be found by the substitution a Tr/2 — 
a.) 




1 .2,0.4) 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
1.0 0.25 0.5 0.75 1.0 1.25 1.5 

alpha 



FIG. 3: Kinematic constraints for strings of unequal ten- 
sions normalized to /i = 3 for the formation of x- and 
z links. The allowed region is to the left of the a::-link 
curves. For z-Iinks, it is above the horizontal line (which 
is absent if ^3 > 2/^2)- The values of /is are again color- 
coded. The values of , 112 , fJ-s) shown are: (1.0,0.6,1.4), 
red, dashed; (1.3,0.3,1.4), red, dotted; (1.4,0.2,1.4), red, 
solid; (1.0,0.8,1.2), blue, dashed; (1.3,0.5,1.2), blue, dot- 
ted; (1.4,0.4,1.2), blue, solid; (1.2,0.8,1.0), green, dashed; 
(1.3,0.7,1.0), green, dotted; (1.4,0.6,1.0), green, solid; 
(1.4,1.2,0.4), magenta, solid. 



When /i_ 7^ 0, it does not seem to be possible to solve 
analytically for all a. We note, however, that there 
arc certain limiting cases for which an analytic solution 
is possible. Firstly, in the low-velocity limit, v ^ 0, it 
is clear from ([55)1 that also u ^ 0. In this limit, sin 6' = 
(p-/ fi^) sin a. It is then straightforward to show that an 
x-link is only possible if 



sm a(„=o) < 



M4 



(45) 



Note that the value of a increases if /is decreases, or if /t_ 
increases. [Recall that the triangle inequalities impose 
the restrictions ((28|) .] 

We can also find a solution for a = 0. Here u is readily 
obtainable from ([55)) and then (j^O)) yields 



(46) 



Thus Vc defined in ()46p increases with fp_. Both these 
limits indicate that in general the kinematic constraints 
exclude a smaller region of the {a,v) plane as the string 
tensions become more different. This behavior is readily 
seen in Fig. [31 
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In the equal tension case, /i_ = 0, Vc is always less 



than 1. Here, however, < I only if 



(47) 



If this condition is satisfied by the tension of the poten- 
tially linking string, there is a velocity Vc above which 
abelian strings will necessarily pass through each other 
rather than intercommuting. However, when the condi- 
tion is violated, the entire high-velocity region v 1 
for all a is included in the allowed region. The bounding 
curve in this case bends to the right, and reaches a = tt/2 
at a finite velocity, above which an x-link is kinemati- 
cally allowed for any angle a. This velocity constraint 
for a = 7r/2 is 



(a=7r/2) 



> 



(48) 



Note that this limit is equal to 



Let us start by looking at the distribution of the vari- 
ables Cj , assuming that the unit vectors are randomly 
distributed on the unit sphere, and mutually indepen- 
dent. Let us choose the z-axis along the direction of bg, 
and bj in the {x,z) plane; h'2 = (sin 61, 0, cos 0). Then 
b'j^ = (sin/3cos 0, sin/3sin 0, cos /?) and we may assume a 
uniform distribution in the variables ci , C2 and cf) where 



C3 = C1C2 + \/ 1 - q W 1 - C2 cos ( 



(49) 



Our aim is to calculate the probability distribution for, 
say, the rate si at which the first string grows. Specifi- 
cally, let P{si)dsi be the probability that si lies between 
si and si + dsi. Clearly, 



/iMi(l-ci) 



X(5 si - 1 + 



(50) 



VI. RATE OF CHANGE OF STRING LENGTHS 

One of the important reasons for studying the kinemat- 
ics of string collisions is that the results may throw some 
light on the question of how a network of such strings 
would evolve in the early universe. If we ignore the Hub- 
ble expansion and any energy loss mechanisms, then the 
energy in the string network is fixed, but some strings 
will shorten and others will grow. We may ask how fast, 
on average, is this growth or shortening. 

It is reasonable to assume that at any string junction, 
the unit vectors b^ representing the ingoing waves are 
randomly distributed on the unit sphere, and mutually 
independent. (This might not be true if for example two 
of the strings come from the same other vertex, but that 
is presumably not a common situation.) If the strings are 
all of the same tension, then because of energy conserva- 
tion it is clear that (sj) must vanish for each but this is 
not necessarily so if the tensions are different. And as we 
shall see, even for the equal-tension case, the zero mean 
does not mean that the distribution is symmetrical; it is 
actually not true that strings are as likely to grow as to 
shrink. 



where we have used ([TOl). and M = M (01,02,03) is de- 
fined in PT|) . 

The resulting calculation is fairly lengthy, so we rele- 
gate the details to Appendix B. The distribution P{si) 
takes three different analytic forms in three intervals, 
namely 



1. -1 < Si < 



9 9 2 

2^1/^2 



2/-tiM3 



2. d^ii^ < < d^ii^. (51) 



(Without loss of generality, we have chosen /i2 > Ms-) To 
write the final expressions for the probability distribution 
concisely, it is useful to introduce the constant 



= 



(52) 



3i^lI^2J^3 

Then the expressions for P{si) in the three regions are 



1. P{h)^Q^i'f^^ 



2 2/l+-^A (2/^2 + 2^3 - /ii) -I- (2^1 - /i2 - /X3)si 



1 - Si 



(/ilSl + ^I2+ Ms)^ 



2. Pis,)^Q.I^^^'+^''-^'^- 



Mi(l-si)2 



(53) 



o „f X ^2 2 /^l + 2/i2 + 2/X3) -f (2pii Al2 + /^3)Sl 

3. P(si) Q^i^/i — — ■ -3 . 

I 



The probability distribution has kinks at the boundaries between these regions, given by (j53|) . The form of the 




FIG. 4: The distribution of P(si) plotted against si. The 
values of fii are indicated by colors. The curves shown are for 
(/ii,^2,M3) = (1.4,0.8,0.8), red, dashed; (1.4,1.2,0.4), red, 
dotted; (1.4,1.4,0.2), red, solid; (1.0,1.0,1.0), black, dotted; 
(1.0, 1.4,0.6), black, solid; (0.2, 1.4, 1.4), green, solid. 



distribution is illustrated for various cases in Fig. |4l 

In the particular case where the tensions are all equal 
the distribution has a single kink, at si = 
is 



i. Its form 



1. P{si) 
3. P(.si) 



27 



1 + Sl 



(si+2)3 VI 
4si + 5 
i-si + 2)3 



Si > 



Si < 



It is interesting that although the mean value of si in this 
distribution is zero, as it must be, the distribution is not 
symmetrical — the dotted black curve in Figure [H At 
any particular time, it is most probable that one of the 
three legs is growing, while the other two are shrinking 
(of course at a slower rate). 

More generally, if fJ-2 ~ fJ-s, there is only a single kink; 
three of the curves in Fig. (|4|) are examples of this case, 
with /ii = 1.4, 1.0, and 0.2. 

It is now straightforward (if tedious) to compute the 
average value of si. We find 



3/ii - /i ^ _Q 
A*i 
z/3)i/|ln 



3/ii 
-(/^iH 



V2)vl In 



-(/ii + 



/i)i^i In 



4//2M3 



4^1^2 



(54) 




FIG. 5: Average value of si plotted against |/i2 — /isj for 
^ti = 0.2 (blue), 0.6 (green), 1.0 (black) and 1.4 (red). 



Note that the expression ([51)1 is symmetrical under the 
interchange ^2 /^3, as it should be. It is also easy to 
check that it satisfies the consistency condition 



Mi(si) + ^2(52) + /U3(s3) = 0. 



(55) 



Finally it is interesting to note that in general (si) is 
more likely to be positive if jii is small, or if the other 
two tensions are very different. In Fig. [5l the average 
value is plotted against |/i2 — /i3| for various values of jjLi. 
Notice that for yUi < 1, (si) is always positive. It appears 
that in a network of strings there may be tendency for 
the lighter strings to grow at the expense of the heavier 
ones. This seems to be consistent with results obtained 
by stud yin g the statistical mechanics of such networks of 
strings [42l |. 



VII. AVERAGE STRING VELOCITY 

For ordinary Nambu-Goto strings in flat space-time, 
the r.m.s. transverse string velocity in a random tangle 
of strings is 1/V2. It is interesting to ask whether this 
figure would be different for junction-forming strings. It 
is easy in principle to answer this question. For example 
the square of the velocity of the first string adjacent to 
the vertex is 



K), 



(56) 



and using dSl) and ^ this can be expressed in terms of 



the Co . We have 



(1 + si)a; • b'l 



(1-ii) 



(1 - S2)C3 



9 



-^^(1 _ S3)c2. (57) Using ([10]) this becomes 
I 



b' 



lyiMijl - ci) - 2/xiM2(l - C2)e3 - 2^1^3(1 - 03)02 
-t/iMi(l - ci) + 2^iAf2(l - C2) + 2^1^/3(1 - C3) 



(58) 



We can then compute the average over the probabihty 
distribution already obtained by plugging this expression 
into an integral of the same form as ([50)1 . In this case, 
there does not seem to be any obvious way of obtaining an 
analytic result, but it is possible to make some progress. 
It seems to be slightly easier to compute the average of 
x'^ rather than x^: 



^^1 / 



2/^2(1 „C2)(1-C3) 



-Mi(l - ci) + 2^i[i/3(l - C2) + Ml - C3)] 

(59) 

Then if we use as independent variables 02,03, -0 with 



ci = C2C3 + ^l-cl^l-clcosi), 
we can perform the ip integration, to yield 

(1-C2)(1-C3) 



(iC3 



(60) 



(61) 



where 

F = 2/xi[i/3(l - C2) + V2{1 - C3)] - V2V^{1 - C2C3), (62) 
and 



G = V2V3\ I 



1 



(63) 



Remarkably ([61]) can be evaluated to give the general 
result 



(x?) = ((l-xi^)) 

13/12 (4^ii -/2)(/ii +/2)2 



i5(Mf - n^) 

^ (4/xi + /2)(mi -/z)^ 2/ii 



15^i(/ii - /z)2 
2 



In 



2Aii 



15/ii(^i + /i)2 



Ml - M 



Ml + M 
(64) 



where we have introduced /2 = /i2 — /X3. A number of 
interesting features can be seen. First, the result is in- 
dependent of ^2 + M3 7 depending only on the single ratio 
Secondly, the limit /l = gives 



1 

15 



[1 + 8 In 2] - 0.436, 



(65) 



0.67 — 



0.66 




FIG. 

1m2 - 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0.25 0.5 0.75 1.0 1.25 



6: R.m.s. value of the string velocity v plotted against 
/ia] for several values of /ii. 



independently of jii. In particular we note that, even for 
strings of equal tension, the r.m.s. velocity is not l/\/2 as 
it is for ordinary Nambu-Goto strings without junctions. 
Surprisingly, on the contrary, this value is obtained when 
/2 = ^1, a limit in which the triangle inequality is only 
just satisfied. 

For completeness, in Fig. [H] we have plotted Wrms = 
{x-i), for various choices of /xi, a function of |m2 — Msl- 



VIII. CONCLUSIONS 

In an earlier letter (ssj . we considered the kinematic 
constraints on the possibility of intercommuting of strings 
that form junctions. Concentrating on the case where 
the approaching strings have equal tensions, we showed 
that if the relative velocity with which they meet is too 
large, no exchange can take place and the strings will pass 
through one another (or, for non-abelian strings, become 
joined by a string in the direction of the relative velocity 
or form a linked X configuration). In this paper we have 
extended the analysis to the more realistic case where 
the approaching strings have different tensions. For some 
values of the string tensions, it may happen that ultra- 
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relativistic strings can exchange partners. 

We first studied the reflection and transmission of 
small-amplitude waves at string junctions — something 
that may be important in predicting the gravitational ra- 
diation from strings. We determined the fraction of en- 
ergy transmitted and reflected across the junction, and 
showed how the reflected wave is generally tilted in the 
opposite direction to the incoming wave. 

In sections iryl andlVl we generalized the kinematic con- 
straints (dH) and ([1^ of I to the case in which the two 
colliding strings have different tensions /xi and fj,2- The 
question of whether strings always intercommute is vital 
in understanding the evolution of a network of strings; in 
particular it affects the final density of strings found in 
the scaling regime of a network, if indeed a scaling solu- 
tion is reached. We have established the criteria required 
for intercommuting in terms of the incoming velocity and 
angle of approach, as a function of the string tensions. A 
particularly important combination of tensions is that 
given in the inequality (j47[) . If this inequality is sat- 
isfied, for instance when fii = /i2, then ultra-relativistic 
strings cannot intercommute. If it is violated, then ultra- 
relativistic strings may intercommute. Our result that 
intercommuting does not always happen appears at first 
sight to be in contrast to the claim in ^l| , but we believe 
the regime where their results are applicable corresponds 
to low velocities of approach as the moduli approxima- 
tion they use breaks down for high velocities. It is in the 
high velocity regime that intercommutation may break 
down. 

Since the initial velocity of colliding strings is crucial 
to determining the importance of the kinematical con- 
straints, in sections IVTl and IVIII we considered the global 
properties of junctions and strings in a network. A plau- 
sible assumption for a string network is that the incoming 
waves at a junction are randomly distributed. This has 
allowed us to make progress in determining the r.m.s. 
velocites of strings. In particular (ignoring energy loss 
mechanisms such as expansion of the universe) we have 
calculated the average speed at which a junction moves 
along each of the three strings from which it is formed. 
Our results are intriguing. For example, even for the case 
of equal tension strings, although the average velocity of 
the junctions is zero as expected, the distribution of the 
velocities is not peaked around zero, but around a nega- 
tive velocity (actually around = —1/2) indicating that 
even in this apparently symmetric case, it is most prob- 
able that at any particular time, one of the three legs is 
growing, while the other two are shrinking, all be it at 
a slower rate. The case of unequal tensions can also be 
solved analytically and our results suggest that junction 
dynamics may be such as to preferentially remove the 
heavy strings from the network. It thus seems likely that 
the system will evolve to one where the lightest strings 
are dominating the dynamics, though of course junctions 
are still present. Regarding the r.m.s. velocities of the 
strings themselves, we showed that they arc generically 
smaller than the 1 /^/2 characteristic of Nambu-Goto net- 



works, even in the case when the strings all have equal 
tensions. 

In a future publication we intend to present some ex- 
act solutions for loops containing junctions. Given the 
new features we have uncovered for the dynamics of 
string networks with junctions, we should not be sur- 
prised to find important results concerning the distribu- 
tion of kinks and cusps in these more complicated con- 
figurations. This in turn could have a bearing on the 
gravitational radiation emitted from loops of strings with 
junctions. 
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APPENDIX A 

Here we outline the calculation of equation ([38|). Let 
/i± = ^1 ± ^2, so that 

Then the second equality in ([37]) yields 
u 

V 

/x_(l + 7^-^7j7-^ cos a cos 6') + ^^^{uv — ^^^^^^ sinaBinO) 
M3(l + l^^lu^ cosacosfl) + fi-{uv — ^~^ju^ sinasinfl) 

Next we can eliminate cos 6 and sinO. ^From the first 
equality in (|57|) . 

^ V cos a 

cosO = , 

V 1)^ cos^ a + sin^ a 

u sin a 

sva.6 = —j^^^^^^^^^^^=. 

V cos^ a + sin^ a. 

Thus multiplying up and grouping all the terms involving 
7"^ on the right, gives 

7-I7-1 „ 

" [{li-v — IJ3u)v cos a + 

V v'^ cos^ a + u'^ sin^ a 

(li-u — fj,3v)usm^ a]. 

Now square to obtain 

(w^ cos^ a + sin^ Q:)[/i3u(l — v'^) — ~ "^)]^ = 

-(1 - v^)(l - u^)[^-(w^ cos^ a + u'^ sin^ a) - Hauv]'^. 
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On expansion, the terms in uv^ which come from the 
cross terms in each square bracket, cancel, and we are 
left with 

F [/^3U^(1 - u^) - /x^(l - u^){v^ cos^ a + v? siv? a)] = 
where 

F = (1 - v^){v^ cos^ a + u^ sin^ a) - ^^(l - u^) 

which vanishes only v — ±u. Thus we finally get the 
quadratic ([38|) for v?: 



li^_{siT? a)u^ + [/X3(l - w^) + ^^(w^ cos^ a - sin^ a)]u^ 

2 2 2 n 

— ^jl_v cos a = 0. 

This always has one positive root as the discriminant is 
positive. 



APPENDIX B 

To carry out the integral in Eq. (jSOp , it is convenient to 
go over to a more symmetrical form by changing variable 
from (j) to C3, using (j49p . This introduces a Jacobian 
factor, which is the inverse of 



(ica 



1 — c{\ I — c\ sin (/) = — 



where 



J = (1-C2)(1-C2)-(C3-C1C2)' 



1-c 



2C1C2C3. 



Note that the physically allowed region in (ci, 02,03) 
space is characterized by J > 0. There is also a fac- 
tor of 2 because there are two values of for each C3. 
Letting tui = 1 — si, we arrive at 



P(l — u)i) = — J dci j dc2 J dc. 



e{j) 



x5\wi 



l-iiM 



(66) 



It is straightforward to perform the C3 integral using 
the delta function. This gives 



F(l — Wl)^- — ^ / rfCill — Ci) / rfC2 



0{J1 

(67) 



where now 

J = (1 - c?)(l - c^) - [c3(w;i) - ciC2]^ (68) 

Here C3 (wi ) is obtained by equating to zero the argument 
of the delta function in (j67p and solving for C3: 

c,[w,) = l + ^(l-ci) (1 ^)+^(l-C2). (69) 
Ma V /ill"!/ ^^3 



Now since, by ([68|) . J is a quadratic function of C2, we can 
go on to perform the C2 integral in (|67p . It is clear from 
([68|l that J < when C2 = ±1, so the effective limits of 
integration for the C2 integral are given by the roots of 
J = 0. If there are no real roots, the contribution is zero. 
Specifically, J has the form 



where 



with 



where 



J ^^{A + 2Bc2-Kcl), 



A={l-cl)Ml-H\ 

B = {M2+ciM3)H, 

K = M| + + 2C1M2M3 



iJ = M2 + M3-(l-ci)r, 



(70) 



Y = Ml 



- 1 



(71) 



(72) 



(73) 



Note that by the triangle inequality, Y is positive. The 
discriminant, which determines whether real roots exist, 
is 



AK={l-cl){K-H') = il-cl)D, 



(74) 



say. If the roots are C2±, then the C2 integral reduces to 



Thus 

P(l - wi) 



(75) 



(76) 



Now from (f7T|) - ((73|) 
D 



= [Ah + Msf - 2(1 - ci)M2M3 

-[AI2 + Ah - (1 - ci)Y]^ 
= (l-ci)[-2Af2M3 

+2{M2 + M3)Y - {1 - ci)Y^]. 



(77) 



It is useful to change variable from ci to ^ defined by 



= K = {AI2 + Ahf - 2(1 - ci)M2A'h 



(78) 



The range of values of ^ corresponding to — 1 < ci < 1 is 
then 



IM2-M3I <i<M2 + A'h. 



(79) 



For convenience, we shall assume in what follows that 
M2 > M3 7 so the lower limit is AI2 — M3. The effect of the 
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condition that the square bracket in ((77)) be positive is 
to require that 

i>ic = ^—^ (80) 



Hence 



Y 



/iMi idi {M2 + M3)^~e 



4.fiiwl M2M3 2M2M3^ 



where 



(81) 



(82) 



C+ = M2 + M3, 
e_ = max(^„Af2- A/3). 

The integral (|5T|) can be rewritten as 



and hence evaluated to yield the final result 



P(l-Wl) 



24/iiM|M. 



2 A/r2 



(84) 



We must distinguish three different ranges of values of 
wi corresponding to the following ranges of Y: 

1. < y < M3, = Y ~ (^"^2 + ^'^3^' 

2. Af3 < y < M2, e- - Af2 - M3; (85) 

2Af2Af3 



3. A/2 < Y, 



(A/2 + A/3) 



Y 



It is straightforward to rewrite the distribution in terms 
of the string tensions /ij and the variable si. The values 
of si corresponding to the three ranges are easily seen to 
be those in (|5ip . and the corresponding expressions for 
P(si) are given in ([53|l . 
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